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1 Summary 

This paper defines a model that could be used to evaluate the ef- 
fectiveness of control measures during the West Africa Ebola epi- 
demic. The introduction gives a general background of the out- 
break, an overview of the disease in context, and common modes of 
regional travel. 

A brief introduction of compartmental models is given, em- 
phasizing parameters useful in modeling. Two major factors are 
addressed: the declining effectiveness of traditional control and the 
possible benefit of increased monitoring and intermittent house- 
level quarantine. 

2 Purpose 

This model has been formulated to aid implementation of the Sierra 
Leone Government's control measures, specifically the develop- 
ment of a national surveillance program combined with rolling 
lockdown. 

Although a large scale Ebola epidemic has never been controlled 
in this manner, the hypothetical effect may be inferred from the 
fact that the host species and humans have historically cohabited 
without a widespread epidemic, or reported outbreak, in West 
Africa. 1 

Further, experience has shown that calls for international assis- 
tance are not necessarily heeded in a timely or effective manner, 
and that the capacity to stop the spread of the disease through con- 
ventional measures is not scaling with the course of the epidemic. 



1 Derek J. Taylor, Matthew J. Ballinger, 
Jack J. Zhan, Laura E. Hanzly, and 
Jeremy A. Bruenn. Evidence that 
ebolaviruses and cuevaviruses have 
been diverging from marburgviruses 
since the miocene. Peer], 2:e556, 
9 2014. ISSN 2167-8359. doi: 
10.7717/peerj.556. URLhttp: 
//dx. doi.org/10.7717/peerj .556 
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3 Background 

The first identified case died on December 6 th , 2013, in Meliandou 
Village, Gueckedou Prefecture, Guinea. 2 Before being identified 2 Baize et al., 2014 

as Zaire Ebola virus (ZEBOV), the outbreak spread to Macenta pre- 
fecture around March 7 th . The disease was identified by Medecins 
sans Frontieres (MSF) from samples taken in mid-March. In early 
April, the outbreak spread to Conakry, a major urban center and 
the capital of Guinea. In early June, outbreaks were reported in 
Telimele and Boffa prefectures in Guinea, Kailahun and Kenema 
Districts in Sierra Leone, and Lofa and Montserrado counties in 
Liberia. 

The outbreak is the largest cluster of Ebola cases on record. Al- 
though previous outbreaks of ZEBOV have seen fatality rates as 
high as 90%, a rough figure of the fatality rate of the this strain of 
ZEBOV is 60%. The incubation period is approximately five days, 
but can range from 3-21 days. The infective period is approximately 
five days. The greatest risk of infection is from very sick individuals 
and the dead. 

3.1 Ebola in West Africa 

Possible reservoirs include fruit bats of the species Hypsignathus 
monstrosus, Epomops franqueti, and Myonycteris torquata, which have 
historically existed in a belt across the whole of West Africa. Ge- 
netic analysis performed on samples from 78 individuals in Sierra 
Leone showed no evidence of additional exposure from animals 
sources, beyond the first case. 3 

Phylogenetic analysis 4 estimates the mutation of this strain oc- 
curring sometime around 2002. Recent analysis of samples collected 
at the Lassa Diagnostic Laboratory in Kenema between 2006-2008 
showed that 8.2% (18/219) of patients received with acute illness, 
having tested negative for malaria, tested positive for EBOV (Zaire- 
76 strain) activity. 5 This data should be qualified by the fact that 
the data is a non-random sample. 6 That is, despite being drawn 
from 219 samples, the data may not be indicative of the general 
population. 

A study of immune responses in epidemic and non epidemic 
regions of Gabon was published in 2010. 7 The study sampled 4,349 
individuals in 220 randomly selected villages, finding a 15.5% av- 
erage of specific anti-ZEBOV IgG immune responses. Spatial varia- 
tion of the resistant fraction correlated strongly with landuse, with 
the highest levels (19.4%) found in forested areas. It was speculated 
that the correlation is the result of conduciveness of the environ- 
ment to bats, and that some non-lethal exposure is transmitted 
from fruit. Specifically the authors note in the discussion: "Thus, 
it is possible that ZEBOV antigenic stimulation or aborted infec- 
tion could occur when villagers handle and eat fruits which have 
been contaminated by bat saliva that may contain infectious virus, 



3 Gire et al., 2014 
* Dudas G, 2014 



5 Schoepp et al, 2014 

6 rather than random, it is individuals 
who sought care at a health clinic, 
showed signs of hemorrhagic fever, 
followed a referral, and then tested 
negative for malaria 

7 Becquart et al., 2010 
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inactivated virus, or simple viral antigens. Brief contact with non- 
infectious viral particles is sufficient to induce the type of specific 
immunity observed in our study." The study also noted a linearly 
increasing immune response with age, indicating that children may 
be particularly vulnerable. The antibodies most commonly found 
in the study were active towards the proteins (VP40 & NP), which 
are multi-functional components critical in reproduction and self- 
assembly of the virus. 



3.2 Travel 

On the available mobility within the region, the road network in 
Guinea is well developed by regional standards, both in terms of 
road condition and overall system implementation. The road net- 
works in Sierra Leone and Liberia suffered from neglect during the 
conflict of the 1990s and lingering set-backs to governance, civil- 
society and investment in the decade following. A number of major 
road infrastructure projects were on-going in Sierra Leone and 
Liberia to improve access to rural areas. The outbreak is currently 
showing the most alarming growth in Liberia, in Lofa County and 
the capital Monrovia. Controlling transmission in a mobile popu- 
lation, located in a remote area, with a high level of cross-border 
trade has proven challenging. 

Public transport is typically provided by modified mini-buses 
carrying up to 21 people in space intended for eight. Fuel subsidies 
make travel accessible to average citizens, with travel across an ur- 
ban area in this manner costing as little as $0.25 USD, and regional 
transport (Freetown-Kailahun) costing around $15 USD. 

In rural areas, inexpensive ($1,250) small displacement Chinese- 
made motorcycles have become increasingly ubiquitous. They can 
carry three to four people, market goods and even livestock. A ten 
kilometer trip costs a single passenger about $1.25 USD one way. 



3.3 Compartmental Models 

Basic compartmental models describing the transmission of com- 
municable diseases as grouped populations were first outlined by 

Kermack-McKendrick 8 . 8 Kermack and McKendrick, 1991 

In brief, compartmental models attempt to quantify the num- 
ber of individuals in various phases of infection with differential 
equations. 

Legrand's Susceptible-Exposed-Infectious-Hospitalized-Funeral- 
Resistant (SEIHFR) model 9 is a compartmental model developed 9 Legrand et al., 2007 

specifically for Ebola. Parameters for the transmission of different 
strains of Ebola virus (EBOV) in different historical situations were 
also developed. Of key interest is the transmission coefficients for 
funerals, community activity and hospitals. Specifically, assuming a 
two day traditional burial period, the contact factor for traditional 
funerals was roughly eight times higher than transmission in the 
community and hospitals combined. And that after effective med- 
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ical intervention, the transmission coefficients for hospitals and 
funerals is effectively zero. 

4 SEIR as Partial Quarantine with Pulse Control 

A simplification of Legrand's model for the current situation would 
be to consider a Susceptible-Exposed-Infected/Funeral-Resistant model, 
where quarantined individuals move to the resistant class but a 
fraction proceeds with an unaltered contact rate. The infective class 
here encompasses both individuals who are infected, as well as 
those who succumb to the disease and are buried in a traditional 
manner. The infection rate approaches a limit as the unquarantined 
infected fraction increases. 

Taking both quarantined and unquarantined cases together re- 
sults in a reproduction ratio above zero but less than the uncon- 
trolled rate, ^controlled ~ 0, ^controlled ^ ^-effective ^ ^uncontrolled- 

Legrand found a basic reproduction ratio of 2.7 (95% CI 1.9-2.8) 
using a stochastic model. However, fitting a deterministic SEIR 
model to the initial outbreak indicates an ^uncontrolled of 1-18, 
with an initial S^ e ff e ctwe or 0-9 after interventions by MSF. This 
may be attributed to the nature of the transmission network and 
be better explained by recent methods developed for hybrid net- 
work /compartmental analysis. 

Following spread to Sierra Leone and Liberia, the & e ff ec tive in- 
creased to around 1.24 overall for several months. After the collapse 
of control measures in Monrovia, the regional overall 2% e fj ec tive has 
increased to ~ 1.8 — 1.9. This model estimates the progression of 
the infection rate as a logistic regression, where /3 approaches some 
value fimax as a function of the infected population. 

Continued transmission implies some unquarantined fraction. 
Unreported cases are not counted in this model. However, the 
unquarantined contact rate together with the current rate of new 
cases (/SSI) may be used to estimate the size of the unreported 
population at some prior time which would lead to a later rate of 
new cases. While, the true extent of unreported cases may be that 

of compounding series 10 , which is not yet treated in this approach. io £~ =Q Ml _ fi_ _ {\ ~S 




In addition to the above assumptions, changes in public behavior 
are difficult to estimate and quantify. And spatial variation can 
have marked effect on contact rates, e.g. between urban and rural 
environments. The progression of the epidemic at different rates in 
different environments presents a complex hazard. 

Using the standard Susceptible-Exposed-Infected-Resistant SEIR 
model, with a separate class for deaths, 
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S' is change in susceptible 
E' is change in exposed 
where: I' is change in infected 
R' is change in resistant 
D' is change in deceased 

and, 



S is susceptible population, 
E is exposed population, 
I is infected population, 
R is resistant population, 
D is deceased population. 



e is the mean exposed period, expressed in inverse days. 
a is the mean infective period, expressed in inverse days. 
/ is the fatality rate. 

/S is the infection rate proportionality constant, which quantifies 
new infectious contacts per infected person per day. With partial 
patient quarantine and tracing measures in place, the effective 
infection rate is assumed to approach the uncontrolled rate as some 
function of I. A binary operator is added for national quarantine 
measures. 



P = (1 - Th) x - 



+ me 



nl 



(2) 



where: 

T is a binary operator for national house-level quarantine. 
h is effect of national house-level quarantine. 
ftmax is the contract factor with public awareness alone. 
m & n are constants. 

and I (as in indigo) is the infected population. 



Begin at the exposure of the index case (late November 2013) 
may be important for fitting the logistic regression equation. A 
simpler or more accurate estimate of the contact rate can certainly 
be substituted for that component. Time driven estimates may be 
fruitful, however the main driving component appears to be the 
number of people infected. 

The desired effect of house-level quarantine and increased con- 
tact tracing is shown in Figure 1. 

A substantial increase in contract tracing should cause a dra- 
matic increase in reported cases, this increase should not to be 
confused with an instantaneous increase in infections (from in- 
dividuals exposed five days prior). The converse should also be 
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Figure 1: Hypothetical effect of inter- 
mittent house-level quarantine on the 
infection rate and the basic reproduction 
ratio, A is the natural frequency of the 
system. 



true. That is, if increased contact-tracing reports fewer cases than 
expected, it would provide evidence that enhancing surveillance 
reduces effectiveness. 

Lockdown measures may increase the contact rate on non- 
quarantine days, due to a higher level of economic activity. The 
system has a natural frequency equal to the exposed period plus 
the infective period. Periods of activity and quarantine should be 
staggered across the natural frequency to avoid erratic behavior in 
the underlying system. A weekly schedule is shown on the timeline 
of Figure 1 . 

Solving for h may show the effect of house-level quarantine 
on the contact rate. If the results are encouraging, a quarantine 
schedule may be developed to minimize the ultimate size of the 
epidemic, as well as the duration. Even if & e ff ec twe does not fall be- 
low 1, mitigating the rate of infection, as soon as possible, increases 
time for development of treatments and vaccines. 

A reliable estimate of fi max should be developed to evaluate the 
time-scale of risks posed by 'wait-for-help' strategies. 
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Acronyms 



ZEBOV Zaire Ebola virus 
EBOV Ebola virus 



2 

3 
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MSF Medecins sans Frontieres 2 

SEIHFR Susceptible-Exposed-Infectious-Hospitalized-Funeral-Resistant . 3 



